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Using Monte Carlo simulations, we compute the integrated emission spectra of black holes in the 
framework of Loop Quantum Gravity (LQG). The black hole emission rates are governed by the 
entropy whose value, in recent holographic loop quantum gravity models, was shown to agree at 
leading order with the Bekenstein-Hawking entropy. Quantum corrections depend on the Barbero- 
Immirzi parameter 7 . Starting with black holes of initial horizon area A ~ 10 2 in Planck units, we 
present the spectra for different values of 7 . Each spectrum clearly decomposes in two distinct parts: 
a continuous background which corresponds to the semi-classical stages of the evaporation and a 
series of discrete peaks which constitutes a signature of the deep quantum structure of the black 
hole. We show that 7 has an effect on both parts that we analyze in details. Finally, we estimate 
the number of black holes and the instrumental resolution required to experimentally distinguish 
between the considered models. 

PACS numbers: 


I. INTRODUCTION 


Loop quantum gravity (LQG) [TH 3 ] proposes a descrip¬ 
tion of the fundamental degrees of freedom responsible 
for the black hole (BH) entropy (see [¥] and references 
therein). According to the most recent results [ 5 H?| , these 
fundamental excitations live on the horizon and are ele¬ 
ments of the Hilbert space of a Chern-Simons theory: the 
gauge group is SU{ 2), the canonical surface is a punc¬ 
tured two-sphere, and the level (coupling constant) k is 
proportional to the horizon area A in Planck units. Punc¬ 
tures are quanta of area associated with horizon-piercing 
edges of the spin-network-states that define the quantum 
states of the bulk exterior geometry. 

As the quantization of a Chern-Simons theory with 
a compact gauge group is now well-understood (see 0- 
m and references therein), the kinematical characteris¬ 
tics of a quantum black hole within the framework of 
LQG are very well-defined. This enables the identifica¬ 
tion of the microstates for a black hole in equilibrium. 
The characterization of their number for a given macro¬ 
scopic horizon area A allows for the computation of BH 
entropy and leads to compatibility with the celebrated 
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Bekenstein-Hawking area law 
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Due to the dependence of the area spectrum on 
the Immirzi parameter 7 compatibility with Bekenstein- 
Hawking BH entropy required fixing the Immirzi param¬ 
eter to a special numerical value in early BH entropy 
calculations [Till [T 2 ]. More precisely in these models the 
BH entropy was shown to be given by 



( 2 ) 


where 70 takes an order one numerical value depending 
on the models mm- Compatibility with |l]) was inter¬ 
preted as a constraint on the value of 7 imposed by the 
existence of the correct semiclassical regime. However, 
the Immirzi parameter is the coupling constant with a 
topological term in the action of gravity mm with no 
impact in the classical equations of motion. Hence, the 
strong dependence of the black holes entropy computa¬ 
tion on 7 remained a controversial aspect. For further 
discussion of this and an exploration of the influence of 
7 in the phase space structure of gravity see EZHBg. 

A promising perspective on the issue of the 7 depen¬ 
dence of BH entropy in LQG was recently put forward 
thanks to the availability of the canonical ensemble for¬ 
mulation of the entropy calculation making use of the 
quasi-local description of black holes [ 20 ]. It was shown 
in m that the semi-classical thermodynamical behavior 
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of BHs can be recovered for all values of 7 if one admits 
the existence of a non trivial ( 7 -dependent) chemical po¬ 
tential conjugate to the number of horizon punctures. 
The problem of the observability of the puncture num¬ 
ber is resolved by taking into account contributions to 
the area degeneracy coming from the matter sector. As 
suggested by the semiclassical description of vacuum fluc¬ 
tuations of non geometric degrees of freedom close to the 
horizon one expects such contributions to degeneracy to 
grow exponentially with the BH area[22]. In [23] one pos¬ 
tulates such a phenomenological contribution and shows 
that compatibility with the existence of large semiclassi¬ 
cal black holes implies that the matter sector saturates 
the holographic bound for a vanishing chemical poten¬ 
tial (the number of punctures drops out, in this way, of 
the list of macroscopic observable quantities). In such 
a scenario, the entropy of large semiclassical black holes 
coincides with 0 to leading order, while the dependence 
on the Immirzi parameter is shifted to sub-leading quan¬ 
tum corrections. 

A possible fundamental explanation the exponential 
degeneracy of the assumption made above was proposed 
in M (and developed further in [25j). In these works, 
the area degeneracy (viewed as an analytic function of 
7 ) is analytically continued from real 7 to complex 7 
and then evaluated at the special complex values 7 = ±i 
to find that it grows asymptotically as exp (A/(41^)) for 
large areas. The result is striking in that these values of 
the Immirzi parameter are special in the connection for¬ 
mulations of gravity: they lead to the simplest covariant 
parametrization of the phase space of general relativity 
in terms of the so-called Ashtekar variables ffil] . These 
results suggest that the quantum theory, when defined in 
terms of self dual variables, might automatically account 
for a holographic degeneracy of the area spectrum of the 
BH horizon. 

Given the present landscape of different models, and 
leaving aside theoretical reasons for preferring one over 
the others, it would be interesting to test their differences 
from an (idealized) observational viewpoint. The value of 
black hole entropy provides little information allowing for 
such a distinction as the models (either by fixing 7 = 70 
or by arguing for an additional exponential degeneracy 
coming from either semiclassical arguments or from an¬ 
alytic continuation to self dual variables) all coincide in 
reproducing ([I]) to leading order. 

The present work shows that a possible way of dis¬ 
tinguishing the models is to analyze out-of-equilibrium 
processes. More precisely, we describe the black hole 
radiation process using statistical mechanical techniques 
together with some assumptions about the dynamical be¬ 
havior of BHs close to equilibrium. In doing this we will 
basically follow the ideas of m leading to the compu¬ 
tation of the integrated emission spectrum for a black 
hole (under some assumptions which will be made precise 
later). The simulated spectrography will allow for the 
study of deviations predicted by the old and new (holo¬ 
graphic) models of LQG from the semiclassical Hawk¬ 


ing process expectation. This is of special interest for 
the holographic models of [231 El] for which the effect of 
7 -dependent sub-leading quantum corrections could be 
detectable. 

The article is organized as follows. In section [TT] we re¬ 
call basic facts about black hole evaporation. In section 
|III| we give a brief description of the different models for 
black holes that we study in this paper. In section |IV| 
we present the numerical methods used for the simula¬ 
tions. Results are presented and discussed in section [V] 
We finish with a discussion in section ED From now on 
we consider units such that G = c = H = 1. In particular, 
areas, lengths and energies are given in Planck units. 


II. EVAPORATION AND TRANSITION RATES 


According to the celebrated no-hair-theorem, station¬ 
ary (axisymmetric) electro-vacuum black holes (those of 
astrophysical relevance) are described by members of the 
Kerr-Newman family and are labelled by their ADM 
mass Af, angular momentum J and electric charge Q. 
The dynamics of the system under small disturbances is 
described by the first law of BH mechanics 

SM = 6 A + n6J + Q6Q, (3) 

07T 

where A is the horizon area, k its surface gravity, $ is the 
horizon electric potential and Q is the black hole electric 
charge. 

Hawking’s semiclassical considerations j2S] imply that 
Kerr-Newman black holes radiate particles according to 
the Planck’s radiation law at temperature T = k/(2tt). 
More precisely, the mean number of particles m 1 q) 

with momentum k a radiated at infinity and carrying en¬ 
ergy w, angular momentum quantum number m and 
charge q is given by 




_ r (u,m,q) _ 

exp (^(w-fim-?$)) ±1 


( 4 ) 


where ui = —k a ^ a and m = k a i/j a are the frequency and 
angular momentum of the mode at infinity, and £ a and 
%l> a are the Kerr-Newman stationarity and axi-symmetric 
Killing vectors respectively normalized at infinity. The 
function T in Q is the so-called greybody factor, and 
the ± sign in the denominator depends on whether one 
is considering radiation of bosons (+) or fermions (—). 

Assuming that when a wave packet carrying energy ui, 
angular momentum m and charge q is emitted the black 
hole parameters are readjusted according to 5M = w, 
SJ = m and SQ = q , we can write, from the first law (|3]) , 
the mean number of emitted particles as follows 


<^> 


r(w,m,g) 
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( 5 ) 


which would make the spectrum Planckian in terms of 
the “area quanta 5A" if it was not for the distortion in¬ 
troduced by T(u>,m,q). This distortion has nothing to 
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do with the thermal properties of the horizon and ap¬ 
pears only when describing the radiation from the point 
of view of observers at infinity: T(ui,m,q) encodes the 
backscattering effects on modes propagating from the 
vicinity of the horizon out to infinity in the background 
Kerr-Newman space-time. 


A. The quasi-local point of view 


From the point of view of stationary observers, denoted 
0 in what follows, at a proper distance l <C M from the 
horizon, one obtains 


<^> 



( 6 ) 


where To is a constant, and thus the spectrum becomes 
strictly Planckian in area quanta. As shown in J5D] the 
(global) first law Q transforms into the following sim¬ 
plest local first law for 0 : 


SE = £ r 6A ’ (7) 

where 8E is the local energy changes as measured by 0. 
The variable Jt denotes the local surface gravity corre¬ 
sponding to the local acceleration of the stationary ob¬ 
servers defined above, i.e. k = ||o ff ||. The (Rindler like) 
near horizon geometry implies that k = 1 /£. 

Notice that k is universal in the sense that it is inde¬ 
pendent of the macroscopic parameters defining the black 
hole [21]. An important consequence of this is that one 
can integrate the local first law to obtain a quasi-local 
energy notion E associated with the (black hole) system 
from the perspective of the observers 0, namely 


E = 


A 

8n£ 


( 8 ) 


Applying the first law 0 to the energy variation corre¬ 
sponds to the emission of a particle with four momentum 
k a and local frequency uj = k a u a as measured by 0 (with 
four velocity u a ) implies that 


8 A 2i tlo 

4 K 


(9) 


which turns ([5]) into a Planckian spectrum at the Unruh 
temperature 


T = 


K 
2l T 


1 

2w£’ 


( 10 ) 


as measured by 0. The frequency independence of To 
in ([6]) in the quasi-local treatment follows from the scale 
invariance of the near horizon approximation 1 . 


1 In the spherically symmetric case, the Klein-Gordon equation 

simplifies to a two dimensional wave equation in the transversal 

coordinates t and r* (where r* is the tortoise coordinate) with 

an effective potential that vanishes at the horizon [29j. This can 
also be shown in the rotating black hole case. 


B. Entropy and transition rates 

The integration of the first law, together with the area 
theorem, leads to the thermodynamical entropy 

S=^+5 0 , (11) 

where Sq is an integration constant that cannot be fixed 
by a thermodynamical reasoning. In fact, as in any ther¬ 
modynamical system, entropy cannot be determined only 
by the use of the first law. It can either be measured 
in an experimental setup (this was the initial way in 
which the concept was introduced) or calculated by us¬ 
ing statistical mechanical methods once a model for the 
fundamental building blocks of the system is available. 
For instance, by computing the microcanonical entropy 
5(A) = log(«yF(A)) where Af(A) is the number of black 
hole microstates compatible with the macroscopic hori¬ 
zon area A. 

When one takes into account the evaporation phe¬ 
nomenon, the black hole can no longer be considered as 
a system at thermodynamical equilibrium. However, for 
dynamical systems which remain close to equilibrium at 
any time during their evolution, the statistical entropy 
is still a well-defined quantity: we say that the system 
is at a local thermodynamical equilibrium. This would 
be the case of large black holes for which Hawking tem¬ 
perature is cold and the horizon area A evolves slowly. 
In such close to equilibrium scenarios, the micro states 
of the horizon compatible with the instantaneous area A 
can be assumed to be equally likely. Thus the probability 
V(A —> A — 8A) for the black hole of horizon area A to 
evaporate a piece of area 8A is well approximated by the 
rate of final to initial number of microstate states 

V (A^A- SA) oc = exp [SS(A)]. (12) 

The proportionality factor is fixed by the requirement 
that the probability density is normalized. Taking into 
account backscattering and normalization issues, one 
shows that the probability for a black hole to emit a par¬ 
ticle whose energy corresponds to the decrease of area 8A 
is given by 

V(A -> A- 8A) = T(A, SA) exp [~SS(A)]. (13) 

We will assume the BH to be of Schwarzschild type. 
Consequently, in the optical limit Mio 1 the greybody 
factors r ss M 2 * ui 2 independently of the particle’s spin. In 
the other limit Mlo <C 1 one has F ss A = 167rM 2 o; 2 for 
scalar particles, T ss 2i tM 2 oj 2 for spin 1/2 particles, and 
r ss yAM 2 u 4 for photons [30]. Numerical limitations 
constrain us to study BHs of a maximum area of about 
10 2 Planck areas. Such BHs have a radius of about 2 
in Planck units which means that the local curvature at 
their horizon is really Planckian. If we boldly stretch the 
assumption of the validity of the semiclassical description 
of space-time outside the black hole, simple dimensional 
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analysis implies that most of the radiated particles fall 
into the Mu <1 regime. 

In order to isolate the quantum effects associated with 
variations of the Barbero-Immirzi parameter it will be 
convenient (in some cases) to assume F(A, 5A) = Fo with 
a constant fixed from normalization issues only. Such a 
choice of T corresponds to neglecting backscattering in 
the description of the emission process and could be re¬ 
garded as describing the physics of the local observers 
0 placed close to the horizon. We will often use this 
terminology in what follows. However, one has to keep 
in mind that while the quasi local perspective is very 
useful in thermodynamical investigations, the notion of 
spectrography does not really make sense for 0. The 
reason for this is that quasi local observers perceive their 
environment at thermal equilibrium at the Unruh tem¬ 
perature and are, in this sense, incapable of discerning 
particles coming from the BH from others in the thermal 
bath. 


III. THREE MODELS OF BLACK HOLES 


According to (13), the evaporation process is com¬ 


pletely governed by the area spectrum of the black hole 
and its entropy. Here we briefly review the main models 
in the literature whose spectrography will be studied in 
section m 


A. Models where 7 must be fixed to 70 

The micro canonical computations taking into account 
only the quantum geometry excitations with no extra de¬ 
generacy associated to non-geometric degrees of freedom 
( e.g. matter fields vacuum fluctuations) lead to the fol¬ 
lowing general expression for the entropy m 

S=^+o(A), (14) 

7 4 

where 70 is a numerical factor of order one depending on 
the details of the state counting [4] . The quantum correc¬ 
tions o{A) are usually logarithmic at the leading order. 
As explained in the introduction, compatibility with the 
first law imposes the constraint 7 = 70 - The model only 
make sense for that special value of the Immirzi param¬ 
eter. There is therefore no interesting physics associated 
with varying 7 . 


B. The holographic models 


matter coupling. If one uses the qualitative behaviour 
of matter degeneracy suggested by standard QFT with 
a cut-off at the vicinity of the horizon (be. exponential 
growth of vacuum entanglement in terms of the BH area), 
then the entropy becomes 


A / 7 T A . 1 — r. 

S= 4 + ^ +0{VI) - 


(15) 


The Immirzi parameter 7 affects only quantum correc¬ 
tions. This result is in agreement with semiclassical ex¬ 
pectations 33| ■ Notice that the entropy has been com¬ 
puted using the quasi-local point of view, the second or¬ 
der corrections to the entropy are explored in [ 321 . 

The presence of the \JA correction can be understood 
directly in the microcanonical framework as follows. In 
the holographic model, the number of black hole mi¬ 
crostates is given by 


JV(A) = exp (-)p(A), 


(16) 


where p(A) denotes the number of ways to write A as the 
finite sum 


A = 8 tt 7 ^2 n iVj(j + !)» 


(17) 


where j are distinct half-integers and n :) are integers. 
At the semi-classical limit (large A), the large spins j 
dominate and then p(A) is nothing but the number of 
partitions of the integer Na = A/ ( 4 ^ 7 ). Its asymptotic 
has been known for almost a century, namely 


logp(A) = 7T 


2N A 

3 


0(\ogN A ) = 



+ 0 (logA), 


which leads immediately to the quantum corrections of 
the entropy |l5|). In the present context one does not 


need to fix the value of the Immirzi parameter and there 
is interesting physics associated to its variations. We will 
study the effects of varying 7 on the BH spectrography 
in section m 

As explained in the introduction, a fundamental un¬ 
derstanding of the holographic hypothesis in the mod¬ 
els leading to (15) might come from the recent results 


3] that indicate a relationship between hologra¬ 
phy and LQG for complex Ashtekar variables (7 = ±i). 
Nonetheless, if in the holographic treatment the degener¬ 
acy originates from the matter fields degrees of freedom, 
it clearly has a (quantum) geometrical origin in the con¬ 
text of complex variables treatment of black holes. One 
can however argue that matter and geometry, which are 
already intimately linked at the classical level, could give 
exactly the same contribution to the black hole degener¬ 
acy (only does the point of view change). 


The previous results correspond to the simplified situ¬ 
ation where purely gravitational degrees of freedom are 
taken into account in the counting. This has been tra¬ 
ditionally the (rather artificial) play-ground for testing 
the theory without worrying about the difficult issue of 


IV. NUMERICAL METHOD 

To study the evaporation process and to analyze the 
effects of the different parameters, we have developed a 
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Monte Carlo simulation procedure. The simulation has 
been adapted from the one used in m- The simula¬ 
tion starts at a given mass, whose value is high enough 
to be well above the deep quantum regime but small 
enough so that all states below this mass can be explic¬ 
itly calculated. In this paper, we start with black holes 
of area Amt = '^ 7r 7min x 70 with 7 m ; n = 0.2 which gives 
Ai n = 1-76 x 10 2 in Planck units, whatever the value of 
7 . This corresponds to an initial BH mass of the order 
of 2 Planck masses, i.e. 


M bh « 2. (18) 

Starting with the same area enables us to make easy 
comparisons. 10 ' black hole evaporations are simulated 
to produce each of the spectra shown in Figure |T| 


In each case, a standard-model particle is randomly 
selected at each step, with a probability weighted by 
its number of internal degrees of freedom and the corre¬ 
sponding (spin-dependent) greybody factor. Only pho¬ 
tons, neutrinos, and charged leptons are finally kept for 
the analysis. Quarks and gluons undergo fragmenta¬ 
tion and generate hadrons with wide energy distribu¬ 
tions. Were new particles beyond the standard mod¬ 
els to be emitted, this would not change drastically the 
results presented here. Only the percentage of mass 
emitted in the photon-neutrinos-charged-leptons chanel 
would change and therefore the number of black holes 
required to reach a given confidence level. But the shape 
of the spectra and the differences between the models 
would qualitatively remain the same. 

Once the particle type is selected, its energy is given by 
the energy loss of the black hole, i.e. E — SM. From the 
(usual) point of view of an observer at infinity, the mass 
M and the horizon area of the black hole are related by 
A = 167 rAf 2 , and then the energy of the emitted particle 
measured by such an observer is given by 

B = = (19) 


where 5A represents the variation of area due to the emis¬ 
sion of the particle. The last identity is true only for small 
variations SA compared to the original area A. The re¬ 
lation between the energy and 5A is simpler from the 
quasi-local point of view. In that scheme, we saw that 
the energy of the black hole is directly proportional to 
the horizon area according to ([ 8 ]). If we denote by Ef 
the quasi-local energy of the emitted particle, with t the 
distance of the quasi-local observer to the horizon, we 
obtain immediately 


Ee 


6A 

8 ttC 


( 20 ) 


From a physical point of view, it is more satisfying to 
plot the spectrum as a function of FA; but it is much 
easier to analyze the results when the spectrum is given 
as a function of Eg. For these reasons, we will give the 


two representations of the black hole spectrum. 

The probability of transition between two black hole 
states of different masses (or areas) is taken to be given 
by the exponential of their entropy difference weighted 
by the greybody factor. Most probably, those greybody 
factors should receive quantum gravity corrections, es¬ 
pecially for the late stages of the evaporation. These 
corrections may become very large and considering only 
a semi-classical greybody factor may then lead to over¬ 
simplified physical predictions. The black hole can un¬ 
dergo a transition to any state having a lower mass. In 
the quantum gravity model we are considering, only dis¬ 
crete values of the area are possible, and the probabil¬ 
ity is driven by the entropy associated with the number 
of states. To make comparisons with the semi-classical 
black hole evaporation, we have also computed numer¬ 
ically the Hawking spectrum which is governed by the 
classical Bekenstein-Hawking area law. 

It should be noticed that, in quantum gravity models, 
the transition to the last state, that is M = 0, naturally 
happens. In the semi-classical Hawking case, however, 
the standard formula has to be slightly modified. The 
Hawking law naively applied would lead the black hole 
to emit, at the last stage, more energy than it has. We 
have therefore used a truncation and considered that 
the energy of the emitted particle is M each time the 
Hawking spectrum would have led to E > M. This only 
affects the very last emission. 

The Monte Carlo is by construction a random process. 
It means that simulating the same process in the same 
conditions would lead to another spectrum, as in real life, 
due to the fundamentally quantum nature of the underly¬ 
ing physics. The spectra presented here are “mean” spec¬ 
tra corresponding to many evaporating black holes ( 10 '). 
As the number of emitted particles in the deep quantum 
gravity regime is quite small (generically less than 10 ) for 
a single black hole, obviously only a statistical analysis 
can lead to a significant conclusion. To remain realis¬ 
tic, we have added an ‘experimental’ uncertainty on the 
reconstructed energy of the detected particles. If the res¬ 
olution was infinite, a single detection would immediately 
allow to distinguish between the Hawking case (that is 
a continuum of states) and the quantum gravity mod¬ 
els (the distinction between quantum gravity proposals 
having the same area eigenvalues but different number of 
micro-states is anyway more subtle). This is however not 
realistic and we have taken into account a kind of reason¬ 
able detector uncertainty. The results presented depend 
on the number of black holes and on the experimental 
resolution (see figure [b]). 


V. RESULTS: BLACK HOLES SPECTRA 


We concentrate our analysis on the holographic model 
where the number of microstates is exactly given by ( 161 . 
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FIG. 1: Spectrum of a holographic black hole for different values of 7 . The spectrum is represented in linear coordinates as a 
function of E a0 in the graph at the top left; it is represented in logarithmic coordinates as a function of Eao in the graph at 
the top right; it is represented in linear coordinates as a function of A A in the graph at the bottom left; it is represented in 
linear coordinates as a function of A A in the graph at the bottom right. 


At the semi-classical limit, the resulting entropy is (15). 


As the greybody factor is fixed to a constant, the proba¬ 
bility for the black hole to emit a particle associated with 
the decrease 8 A of the horizon area is totally governed by 
the entropy according to (12). In the quasi-local point of 


view ( 201 , such a probability amplitude can be expressed 


in terms of 5A and is directly related to the black hole 
spectrum because the energy of the particle Eg is propor¬ 
tional to 8A. From the point of view of an observer at 
infinity, the spectrum stems from the probability "P(A, E ) 
which, as we will argue below, takes the form 


V{A,E) oc E 2 exp[—(55(A)], 


( 21 ) 


with the constant fixed by normalization issues. 

Using the Monte-Carlo simulation we have just de¬ 
scribed in the previous sections, we obtained: 


1 . exact spectra with different values of 7 of order 
10 -1 (Fig. 0 We considered the following set 
{0.2; 0.3; 0.4} of values for 7 . We have ploted the 
spectra in terms of 8A (or equivalently in terms 
of the local energy Eg) and also in terms of E^. 
We have presented the results both in linear and 
logarithmic representations. 


2 . an accurate comparison of the spectra with the clas¬ 


sical Hawking spectrum. 

3. the smeared spectra which take into account an in¬ 
strumental resolution (Fig[5j). 

Note that the choice of the values of 7 are motivated by 
the value for 7 obtained from old models (7 ~ 0.2). Such 
a choice allows for an easier comparisons with previous 
results. 

Let us start with some general observations. Figure 
[T] shows a strong dependence of the spectrum on the 
Barbero-Immirzi parameter even though the leading 
order term of the black hole entropy (which mainly 
governs the evaporation process) is independent of 
7 . This dependency, which is completely quantum 
gravitational in nature, originates from two aspects: 7 
enters in the discretization of the area spectrum and 7 
shows off in the sub-leading corrections to the entropy. 
Before going deeper into the explanation of the role of 
7 , let us comment further on the general aspects of the 
spectrum. 

The spectrum of emission, as seen from infinity, sepa¬ 
rates into two rather distinct parts: a continuous back¬ 
ground whose amplitude is most important in the in¬ 
frared and a series of discrete peaks which go from the 
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infrared to the deep ultraviolet regime. This separation 
is apparent in the top two panels in figure [l] while it is 
hidden in the bottom panels because the range of the 
local emission plots in terms of AA has been chosen to 
illustrate the continuum region of the spectrum only. No¬ 
tice that the greybody contribution, the factor E 2 in the 
probability amplitude, is responsible for the peak in the 
emission spectrum as seen from infinity in the continuous 
lower energy region. As the greybody factors are trivial 
in the local framework we see a monotonic linear behav¬ 
ior in the logarithmic plot on the bottom right panel. 
Coming back to the bottom panels, we observe that the 
quantum peaks are mainly present in the deep ultravio¬ 
let regime. While the continuous part of the spectrum 
corresponds to particles emitted at the early stage of the 
evaporation process, the peaks reveal the deep quantum 
structure of the black hole corresponding to the latest 
stages of the evaporation when the discrete structure of 
the area is no longer negligible. These two regions of the 
spectrum exhibit structures that are dependent on the 
Barbero-Immirzi parameter. We will now describe these 
features in more detail. 


A. The continuous background: link to the 
semi-classical evaporation 

In both the quasi-local and the infinity points of view, 
spectra have a continuum background which is exponen¬ 
tially decaying. 


1. The quasi-local point of view: analysis of Sp(AA) 

This aspect is best observed in the logarithmic scale 
when the spectrum is viewed as a function of A A. In 
the infrared, the spectrum can be modeled, in a first 
approximation, as follows: 

Sp(AA) ~ Cexp(— tAA), (22) 

where r is the exponential decay rate and C is indepen¬ 
dent of A A. It is interesting to observe that r depends 
on the Barbero-Immirzi parameter (and probably also 
on the normalization constant C but this aspect is less 
interesting). We clearly see that r decreases with 7 . 

All these observations are in fact easy to interpret. In¬ 
deed, in the infrared, the spectrum can be reasonably 
approximated by its semi-classical expression which is 
given by the probability rate V(A —► A — AA). This ap¬ 
proximation immediately implies the following estimate 
for r: 


5S 1 / / 2 tt \ 

T ~ ~5A ~ \ 37 A ) ’ 


(23) 


where we used the semi-classical expression for the en¬ 
tropy (15) up to the first quantum correction. Here A 


can be identified with the initial black hole area, i.e. 


A = 1.76 x 10 2 , as the continuous spectrum concerns 
the early stage of the evaporation process when the black 
hole area remains close to its initial value. The previous 
expression for r is consistent with the numerical simula¬ 
tions. 


We can go further in the analysis of Sp(AA). We notice 
that the exponential decay approximation is clearly bet¬ 
ter in the infrared than in the ultraviolet where, even if 
the statistical fluctuations are more important, the slope 
r seems to increase with A A: the spectrum is steeper for 
large AA (larger than 30 in Planck units) compared to 
small A A (between 0 and 10 Planck units). The previ¬ 
ous expression for r ( [23] ) allows us to easily interpret this 
fact. In the ultraviolet, the area A in (23) can indeed no 
longer be identified with the initial area because we are 
far from the earliest stages of the evaporation process. In 
that case, A has to be identified with the instantaneous 
area of the black hole which is obviously smaller than the 
initial one. The consequence is that the slope r increases 
because it scales as 1/ y/A. The more we go to the ultra¬ 
violet, the larger the slope is, exactly as observed in the 
numerical simulations. 


Let us finish the analysis of Sp(AA) with some remarks 
concerning the comparison with the Hawking spectrum. 
We observe that the quantum and the Hawking spectra 
are closer one to the other in the infrared than in the ul¬ 
traviolet. The reason is simple and comes from the fact 
that quantum corrections to the entropy are no longer 
negligible in the late stages of the evaporation. From the 
expression of the entropy, we see that the larger 7 is, the 
smaller quantum corrections are and then the better the 
approximation to the Hawking spectrum in the ultravio¬ 
let is. This theoretical prediction is consistent with the 
numerical calculations (see Figure [l]). 

In the infrared, the situation is somehow the reverse. 
In that case, the quantum corrections to the BH entropy 
are negliguible, and thus cannot explain the differences 
between the spectra with different 7 . We observe that the 
smaller 7 the better the approximation to the Hawking 
spectrum. This simply comes from the value of the area 
gap in the area spectrum of the black hole which is linear 
in 7 . When 7 is small, the area spectrum of the black 
hole resembles a continuous spectrum and then we expect 
the emission spectrum to be very close to the Hawking 
one, as shown by the numerical simulation. 


2. The point of view at infinity: analysis of Sp(E) 


When the spectrum Sp(-E) is viewed as a function 
of E (which is relevant for the point of view of an ob¬ 
server at infinity), we observe that its continuous back¬ 
ground admits a maximum and the smallest energies are 
suppressed. From a theoretical point of view, this is 
explained because, near the maximum, the continuous 
background can be approximated by its semi-classical ex- 
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pression which leads to the following estimate 2 


Sp (E) ~ CE 2 * exp (-f £), (25) 


where the slope f is now given by 


5S f—r I 87t 2 

T~ — ~2s/^A+ s j—. (26) 


As in the previous analysis, C can be approximated by a 
constant and A can be reasonably identified with the ini¬ 
tial area. From this result, we can immediately estimate 
the energy E max = 2/f for which the spectrum is maxi¬ 
mum together with the value of this maximum which is 
given by 


Sp (E max ) ~ C . (27) 


Therefore, we can predict that the maximum of the spec¬ 
trum increases with 7 . This is exactly what we observe in 
the numerical experiments. Such an analysis shows that 
the quantum corrections to the entropy manifest them¬ 
selves also in the infrared part of the spectrum, and not 
only in the deep ultraviolet regime, what one would have 
expected. This aspect seems to us very interesting to 
underline. 


2 As the initial mass of the BH R 2 2 the continuous part of the 
spectrum falls in the regime Mu <C 1 (putting all the units in 
the previous criterion gives (M (hu /c 2 ))/M 2 = {ME)/M p < 1). 
In the Figure |V A 2| below we show the exact distribution of Mu 
for 7 = 0.4. This means that we can use the low energy form of 
the form factors m from which we would get 

(■ NfdfE 2 + N ph a ph E 4 ) exp (-f E) « 

NfdfE 2 exp (— ?E) , (24) 

where Nf and N p h are the number of fermions and photon de¬ 
grees of freedom and af and 7 ,/, are parameters of the same 
order. The photon term is subdominant as Nf > N p h in the 
standard model. 



FIG. 2: Distribution of Mu for 7 = 0.4 illustrating the fact 
that the continuous part of the spectrum is radiated in the 
IR regime Mu < 1. 


B. The quantum peaks: quantum nature of 
space-time 

In Figure [TJ the linear plot (Sp(AA) or Sp(2?)) is the 
more convenient to analyze and interpret the properties 
of the quantum peaks in the black hole spectrum. These 
peaks (above all the ones with a higher energy) are ex¬ 
pected to appear at the late stages of the black hole evap¬ 
oration process and therefore should correspond to the 
transitions between low area states. In that respect, they 
provide us with a concrete signature of the deep quantum 
structure of the black hole. 

More explicitly, each peak is produced by a transition 
from a state i (of quantum area Ai) to a state j < i (of 
quantum area Aj). Whether we consider the local point 
of view or the point of view at infinity, the energy of the 
corresponding emitted particle is given by 

E Uj = oc 7 , (28) 

Eooij = — \/4j) ^ Vfo (29) 

and therefore it scales respectively as 7 or ^ 7 . This scal¬ 
ing can be immediately observed in the linear plot Sp(<5A) 
and Sp [E). For instance, the three highest peaks in Sp(5) 
correspond to the following approximate energies: 


7 

0.2 

0.3 

0.4 

fo4(7) 

1.60 

2.40 

3.20 

5 A /7 

0.80 

0.80 

0.80 


Indeed, A( 7 ) = 47 r 7 (2y/l x 3 — \/2 x 4), which is the 
area difference between any configuration containing 
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(at least) two j = \ punctures and the configuration 
obtained by replacing them by one j = 1 puncture. 
All other area-difference peaks can be identified to a 
puncture change or removal in the same manner. 

The same exercise for highest energy peaks can be done 
and the observation of the scaling of the energy as y/j 
in the point of view at infinity can be checked. Nev¬ 
ertheless, since SE depends on 5A and A , peaks are 
to be associated with transition between two configu¬ 
rations. For instance, the highest peaks correspond to 
the transition (n,j = ( 2 , 0 , 0 ...) —> (nj) = ( 0 , 0 ,...) 
with Eao = v / 73 1 / 4 2 -1 / 2 , which gives 0.42,0.51,0.59 for 
7 = 0.2,0.3, 0.4 respectively. 

A more delicate effect of 7 concerns the amplitude of 
the peaks. This effect is easier and more precise to study 
in the quasi-local point of view where the amplitudes 
are given by Sp (Euj). In fact, to take into account the 
normalization issues, it is more judicious to study relative 
amplitudes of the form Sp(Euj) /Sp(£/o) where Eq is the 
reference energy corresponding to a given transition iq to 
jo- One can for instance choose the highest energy peak 
as the reference. Such relative amplitudes depends on 7 
according to 

Sp^o) a ex P[ _27r ^(^b ~ £«))]• (31) 

As the energy scales as 7 in the quasi-local point of 
view, the relative amplitude is more suppressed for larger 
7 . This explains the observation that when 7 increases, 
lower energy peaks have larger amplitudes. In the infi¬ 
nite viewpoint, although it is not a priori evident, we 
observe that the presence of greybody factors do not al¬ 
ter the qualitative picture in the range of 7 € [0.2,0.4] 
considered here. We leave the issue of how these features 
extrapolate to a larger range of 7 for further studies, 
the main message here being that the effect of 7 on the 
spectrum (in both its semi-classical and deep quantum 
regimes) is more than a simple dilation. 


C. Comparison with older models 

For completeness we here show in figures [3] and [3] the 
7 dependence of the integrated spectrum of the non- 
holographic models considered in previous literature Ill- 
Notice that a comparison of different values of 7 is only a 
mathematical exercise as in the framework of such mod¬ 
els the Immirzi parameter must be fixed to the value that 
is compatible with Sbh = A /4 which in this case corre¬ 
sponds to 7 « 0.274 [¥]. 

A question of physical relevance is how many evapo¬ 
rating black holes would be necessary to distinguish be¬ 
tween the different models studied here and the Hawking 
semiclassical model. In order to do this in a quite re¬ 
alistic way, we introduced an uncertainty in the energy 
measurement of emitted particles through a random er¬ 
ror. With this modification, the spectra in figure |T| gets 



FIG. 3: 7 dependence of the integrated spectrum in the stan¬ 
dard LQG model (linear scale). 



FIG. 4: 7 dependence of the integrated spectrum in the stan¬ 
dard LQG model (logarithmic scale). 


smeared as shown in figure [5} With such smeared spectra 
is now possible to compute how many detections would 
be necessary to distinguish the different models. Results 
are shown in figure [ 6 ] 

In each case, if the number of black holes is high 
enough and/or the uncertainty small enough, it is possi¬ 
ble, through a Kolmogorov-Smirnov test, to statistically 
distinguish between models. Due to the structure and 
amplitude of the peaks, it is slightly easier to distinguish 
between the pure Hawking spectrum and the old LGQ 
model than between the pure Hawking spectrum and the 
holographic one. 


VI. CONCLUSIONS 

The study presented in this article investigates the 
precise emission from light quantum black holes. Start¬ 
ing with an initial horizon area A ~ 10 2 in Planck units, 
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FIG. 5: 7 dependence of the integrated spectrum in the holo¬ 
graphic model with a simulated error of energy detection: a 
particle of energy E has detected energy E + e, where e is a 
normal distribution with variance a = 0.05 E. 
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FIG. 6 : Number of observed events necessary to distinguish 
the semiclassical Hawking model from LQG models, as a func¬ 
tion of the relative uncertainty of the apparatus. 

the spectra for different values of 7 were computed, 
averaging over many realizations. The continuous 


background corresponding to the semi-classical stages 
is complemented by discrete peaks associated with the 
deep quantum regime. We have shown that the Barbero- 
Immirzi parameter has an important, and sometimes 
subtle, effect on both parts of the spectrum. Finally we 
have calculated the number of black holes, for a given 
experimental uncertainty, required to experimentally 
distinguish between the holographic model and the 
semi-classical Hawking spectrum. 

The analysis presented in this article is more relevant 
at the conceptual level than at the strictly experimental 
level. It is more a “thought experiment”, useful to 
understand subtleties that the analytical investigations 
do not reveal easily, than a proposal for a real experi¬ 
ment. Is it however conceivable that evaporating black 
holes can really be used to probe the models ? In the 
framework of a cosmological production of primordial 
black holes, this is extremely unlikely (see [ 33 ] for a 
review). A very wide mass spectrum of primordial black 
holes is now disfavoured as it would require either a 
high normalization of the primordial fluctuation power 
spectrum (to produce a density contrast of order unity) 
or a blue tilted spectrum, both being excluded by CMB 
observations. There are many other means for creating 
primordial black holes with a narrower mass spectrum. 
In any case, however, the emission spectrum roughly de¬ 
creases with the energy as E~ a with a ranging between 
around 1 below 100 MeV to around 3 above this scale 
[351 SSj. Unless a very strong boost phenomenon occurs 
at late times, it is therefore basically impossible to see 
the last stages without having first detected, at a much 
higher level, the semi-classical stage. A more promising 
avenue would be associated with the production at col¬ 
liders. This however requires a low Planck scale so that 
the collision is trans-planckian and creates a black hole. 
This is not, in principle, incompatible with LQG but 
relies on strong assumptions -such as extra-dimensions- 
that the theory precisely allows to avoid. At this stage, 
the observational detection is therefore unlikely and 
simulations are mostly used as gedankenexperiments to 
understand better the detailed features of the model. 
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